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EXPLICIT RATES OF EXPONENTIAL CONVERGENCE 
FOR REFLECTED JUMP-DIFFUSIONS ON THE HALF-LINE 

ANDREY SARANTSEV 


Abstract. Consider a reflected jump-diffusion on the positive half-line. Assume it is stochas¬ 
tically ordered. We apply the theory of Lyapunov functions and find explicit estimates for the 
rate of exponential convergence to the stationary distribution, as time goes to infinity. This 
continues the work of Lund, Meyn and Tweedie (1996). We apply these results to systems of 
two competing Levy particles with rank-dependent dynamics. 


1. Introduction 

1.1. Exponential convergence of Markov processes. On a state space X, consider a 
Markov process X = {X{t),t > 0) with generator At and transition kernel P^{x,-). Exis¬ 
tence and uniqueness of a stationary distribution vr and convergence X(t) —)■ tt as t —)• oo 
have been extensively studied. One common method to prove an exponential rate of conver¬ 
gence to the stationary distribution vr is to construct a Lyapunov function: that is, a function 
U : X —>■ [1, oo), for which 

MV{x) < —kV{x) + blE{x), X G X, 

where b,k > 0 are constants, and E is a “small” set. There is a precise term small set in this 
theory. In this article, X = M+ := [0, cxo), and for our purposes we can assume E is a compact 
set. If there exists a Lyapunov function V, then (under some additional technical assumptions: 
irreducibility and aperiodicity), there exists a unique stationary distribution tt, and for every 
X G X, the transition probability measure E*(x, •) converges to vr in total variation as t —)■ oo. 
Moreover, the convergence is exponentially fast. More precisely, suppose H-H denotes the total 
variation norm or a similar norm for signed measures on X. (In Section 3, we speciy the exact 
norm which we are using.) Then for some positive constants C(x) and x, we have: 

(1) ||E‘(x,-)-vr(-)||<U(x)e-h 

Results along these lines can be found in III11251ES], as well as in many other articles. Similar 
results are known for discrete-time Markov chains; the reader can consult the classic book [21] . 
However, to estimate the constant x is a much harder task: See, for example, [31 UnilSTl 12811291 
|3T1[39]. In the general case, the exact value of x depends in a complicated way on the constants 
b and A:, on the set E, and on the transition kernel E*(x, •). 

Under some conditions, however, we can simply take x = k. This happens when X = M+, 
E = {0}, and the process X is stochastically ordered. The latter means that if we start two 
copies X' and X” of this process from initial conditions x' < x", then we can couple them so that 
a.s. for all t > 0, we have: X'{t) < X"{t). This remarkable result was proved in [2TI Theorem 
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2.2], (A preceding paper [22] contains similar results for stochastically ordered discrete-time 
Markov chains.) In addition, in [211 Theorem 2.4], they also prove that even for a possibly 
non-stochastically ordered Markov process on M_|_, if it is stochastically dominated by another 
stochastically ordered Markov process on M+ with a Lyapunov function with E = {0}, then the 
original process converges with exponential rate x = k. Let us also mention a paper [SU], which 
generalizes this method for stochastically ordered Markov processes when E ^ {0} (however, 
the results there are not nearly as simple as x = k). 

1.2. Our results. In this paper, we improve upon these results. First, in Theorem 14.11 we 
prove that x = k for stochastically ordered processes (a version of [211 Theorem 2.2]) under 
slightly different assumptions, with an improved constant C{x). Second, in Theorem 15.21 we 
prove a stronger version of m Theorem 2.4] for non-stochastically ordered processes (because 
the norm in ([T]) is stronger in our paper). In particular, our result allows for convergence of 
moments, which does not follow from [211 Theorem 2.4]. Third, in Lemma 16.11 we show that in 
a certain case, this rate x of convergence is exact: one cannot improve the value of x; this serves 
as a counterpart of [211 Theorem 2.3]. Next, we apply this theory to reflected jump-diffusions 
on M+. 


1.3. Reflected jump-diffusions. A reflected jump-diffusion process Z = {Z{t),t > 0) on the 
positive half-line M+ is a process that can be described as follows: As long as it is away from 
zero, it behaves as a diffusion process with drift coefficient g{-) and diffusion coefficient o'^(-). 
When it htis zero, it is reflected back into the positive half-line. It can also make jumps: Take a 
family {i'x)x>o of Borel measures on M+. If this process is now at point x G M+, it can jump with 
intensity r[x) = z/a;(M+), and the destination of this jump is itself a random point, distributed 
according to the probability measure r“^(a;)z/a,(-). If r{x) = 0, then this process cannot jump 
from X. A rigorous dehnition is given in Section 2. 

We can similarly dehne reflected diffusions and jump-diffusions in a domain D C W^. These 
processes have many applications: Among a multitude of existing articles and monographs, let 
us mention [smsi 12D1 EH ESI and references therein for applications to stochastic networks. 
Existence and uniqueness of a stationary distribution and convergence to this stationary distri¬ 
bution as t —>■ cxD for these processes were intensively studied recently. Among many references, 
we point out the articles [211711I2] for reflected diffusions and [TH] [H [23] for reflected jump- 
diffusions. However, these papers do not include explicit estimates of the exponential rate of 
convergence. 

In this paper, we hrst prove a general exponential convergence result for a reflected jump- 
diffusion on M+: this is Theorem 13.21 which does not provide an explicit rate x of exponential 
convergence. Next, we And an explicit rate of convergence for a stochastically ordered reflected 
jump-diffusion in Theorem 14.31 and for a non-stochastically ordered reflected jump-diffusion 
(dominated by another stochastically ordered reflected jump-diffusion) in Corollary 15.31 

1.4. Systems of competing Levy particles. Finally, we apply our results to systems of two 
competing Levy particles, which continues the research from [lElES]. In these systems, each 
particle is a one-dimensional Levy process. Its drift and diffusion coefficients and the jump 
measure depend on the current rank of the particle relative to other particles. Such systems are 
applied in mathematical hnance in [HI H] [16] . 

1.5. Organization of the paper. In Section 2, we introduce all necessary dehnitions, and con¬ 
struct there reflected jump-diffusion processes. In Section 3, we prove exponential convergence 
under some fairly general conditions, but without finding or estimating a rate of exponential 
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convergence. In Section 4, we prove x = A; for stochastically ordered processes, and in Section 
5, for processes dominated by a stochastically ordered process. In Section 6, we show that in a 
certain particular case, our estimate of the rate of convergence is exact. Then we apply these 
results in Section 7 to systems of two competing Levy particles. 

1.6. Notation. Weak convergence of measures or random variables is denoted by We denote 
M+ := [0, cxd) and M_ := (—cxd, 0]. The Dirac point mass at the point x is denoted by Take 
a Borel (signed) measure z/ on M+. For a function / : M+ —)■ M, denote (z/, /) := /dz/. For a 

function / : M_|_ [1, oo), we dehne the following norm: 

(2) ||z/||/:= sup |(z/,^)|. 

l9l<f 

For / = 1, this is the total variation norm: 11-11/= |H|tv- In the rest of the article, we operate on 
a hltered probability space (D, (J^t)t>o,P) with the hltration satisfying the usual conditions. 

For a function / : M+ —)■ M, we let ||/|| := supj,>o |/(a^)|- We denote by Exp(A) the exponential 
distribution on the positive half-line with mean (and rate A). 

2. Definition and Construction of Reflected Jump-Diffusions 

First, let us dehne a rehected diffusion on M+ without jumps. Take functions g : M+ M 
and a : M+ —)■ (0, cx)). Consider an (J^Jt>o-Brownian motion W = (W(t),t > 0). 

Definition 1. A continuous adapted M 4 .-valued process Y = (Y{t),t > 0) is called a reflected 
diffusion on M+ with drift coefficient gf) and diffusion coefficient cr^(-), starting from yo G M+, 
if there exists another real-valued continuous nondecreasing adapted process I = (t(t),t > 0) 
with 1(0) = 0, which can increase only when Y = 0, such that for f > 0 we have: 

Y{t)=yo+[ 9 {Y{s))ds+ [ (T(y(s))dlF(s)-M(t). 

Jo Jo 

Assumption 1. The functions g and a are Lipschitz continuous: for some constant Cl > 0, 

Igix) - giy)\ + flix) - a{y)\ < Cl\x - y\, for all x,y eR+. 

Moreover, the function a is bounded away from zero: inf^-^o c’'(t) > 0. 

It is well known (see, for example, the classic paper [321) that under Assumption 1, for every 
yo G M+, there exists a weak version of the rehected dihusion from Dehnition [T], starting from 
yo, which is unique in law. Moreover, for diherent starting points yo G M+, these processes form 
a Feller continuous strong Markov family. We can dehne the transition semigroup P* which acts 
on functions: / i—)■ P*/, as well as the transition kernel P*(x, •) and the generator A: 

(3) Af{x) = g{x)f{x)+ ^a^{x)f"{x), if /'(O) = 0. 

Take a family {i>x)x>o of hnite Borel measures z/^, on M+. 

Assumption 2. The family {i'x)x>o is weakly continuous: Uxo for Xn Xo- In addition, 

the function r(x) := z/a;(R+) is bounded on M+: sup 2 ,>or(a;) =: p < oo. 

Take a rehected dihusion Y = {Y{t),t > 0) on M+ with drift coefficient g and dihusion 
coefficient Using this process Y, let us construct a weak version of the rehected jump- 
dihusion Z = {Z{t),t > 0) with the same drift and dihusion coefficients and with the family 
(i^x)xGK+ of jump measures, starting from yo G M+. One way to do this is piecing out. 
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For every y G M+, take infinitely many i.i.d. copies n = 1, 2,... of the reflected diffusion 

y, starting from = y. For every x G M+, generate infinitely many i.i.d. copies of 

a random variable ~ r“^(a;)z/a.(-), independent of each other and of the copies of the reflected 
diffusion Y. We assume all processes adapted to (Ji)t>o, and every cx-subalgebra Tt 

contains all for x G M+ and n = 1,2,... Start a process We kill it with intensity 

r(F’^^°’^)(t)): If Cl is the killing time, then 

(4) P (Cl >t) = exp (^- {Y^y°’^\s)) ds 

If Cl < C) 0 , let Xi := F*^^°’^)(Ci), and let yi := Start the process 

C 2 with intensity r{Y^y^’‘^\t)), similarly to (0]), etc. Because r{x) < p ioi x E M+, one can hnd 
a sequence of i.i.d. pi,P 2 , • • • ~ Exp(p) such that a.s. for all k we have: Cfc ^ Vk- Therefore, 

(5) Tk := Cl + ■ ■ ■ + Ck > Vi + ■ ■ ■ + Vk ^ 00 a.s. as k ^ 00 . 

Dehne the process Z = {Z{t),t > 0) as follows: for t G [rfc,rfc+i), k = 0,1,2,..., let Z{t) = 
— Tfc). In other words, it jumps at moments ri,r 2 ,..., and behaves as a reflected 
diffusion without jumps on R+ on each interval (rfc,rfc+i). Because of ([5]), this dehnes Z(t) for 
all t > 0. The following result is proved in [341 Theorem 2.4, Theorem 5.3, Example 1]. 

Proposition 2.1. Under Assumptions 1 and 2, the construction above yields a Feller continuous 
strong Markov process on M+ with generator 

(6) C = A + Af, 

where the operator A is given by ([3|), and Af is defined by 

roo 

(7) Affix) := [fiy) - fix)]u^idy). 

Jo 

3. Lyapunov Functions and Exponential Convergence 

In this section, we dehne Lyapunov functions of a Markov process on M+ and relate them to 
convergence of this Markov process to its stationary distribution with an exponential rate. We 
apply this theory to the case of rehected jump-diffusions. However, we do not hnd an explicit 
rate x of exponential convergence: this requires stochastic ordering, which is done in Section 4. 
Our dehnitions are taken from |33] and are slightly diherent from the usual dehnitions in the 
classic articles [HI |25l |26] . These adjusted dehnitions seem to be more suitable for our purposes. 

3.1. Notation and definitions. Take a Feller continuous strong Markov family (P*)t>o on R+ 
with generator At, which has a domain P(At). Let P*(x, •) be the corresponding transition 
kernel. Slightly abusing the terminology, we will sometimes speak interchangeably about the 
Markov process X = (X(t), t > 0) or the Markov kernel {P^)t>o- We use the standard Markovian 
notation: /iP* is the result of the action of P* on a measure /i; symbols P^, and E^, correspond 
to the copy of X starting from X (0) = x. 

Definition 2. Take a continuous function V : R+ -E [l,C)o) in the domain P(A1). Assume 
there exist constants b,k, z > 0 such that 

(8) JYlVix) < —kV{x) + 61[o,q(a;), x G R+. 

Then V is called a Lyapunov function for this Markov family with Lyapunov constant k. 
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Let us now define the concept of exponential convergence to the stationary distribution. Take 
a function W : M+ [l,C)o) and a constant x > 0. Recall the definition of the norm H-Hpv 
from (|2]). 

Definition 3. This process is called W-uniformly ergodic with an exponential rate of conver¬ 
gence X if there exists a unique stationary distribution tt, and for some constant D, we have: 

(9) \\P^{x,-) — 'n{-)\\w < DW{x)e~'^, a: G M+, t > 0. 


Finally, let us introduce a technical property of this Markov family, which allows us to link 
Lyapunov functions from Dehnition |2] with exponential convergence from Definition |3l 

Definition 4. The Markov process is called totally irreducible if for every t > 0, x G M+, and a 
subset A C M_,_ of positive Lebesgue measure, we have: P*(x, A) > 0. 

The following result is the connection between Lyapunov functions and exponential conver¬ 
gence. It was proved in |33] and is slightly different from classic results of [25l [26l HI] . 

Proposition 3.1. Assume the Markov process is totally irreducible with a Lyapunov function V. 
Then the process is V-uniformly ergodic, and the stationary distribution tt satisfies (tt, R) < oo. 


3.2. Main results. Let us actually construct a Lyapunov function for the reflected jump- 
diffusion process Z = {Z{t),t > 0) from Section 2. We would like to take the following function: 

(10) Vx { x ) := x G M+, 

for some A > 0. Indeed, the first and second derivative operators from (I3|), included in the 
generator C from (jH]), act on this function in a simple way. However, R((0) ^ 0, which contra¬ 
dicts (I3|). Therefore, we cannot simply take Vx as a Lyapunov function; we need to modify it. 
Fix S 2 > Si > 0 and take a nondecreasing function ip : R+ M+ such that 

(11) (p(s) = i ipfs) < s for s > 0. 

[s, s > Sa; 

The easiest way to construct this function is as follows. Take a molliher: a nonnegative C°° 
function a; : M —?• M+ with /jga;(x)dx = 1, with support [—e,e\, where £ = (sa — si)/3. One 
example of this is 

u:e{x) = cexp (—(^ — |3^|)~^) , c > 0. 

Apply this molliher to the following piecewise linear function: 

{ 0, s < Si 

(s2-£)^f5f^, si+e<s<sa-£; 
s, s> S 2 - e. 

The convolution of ip and gives us this necessary function tp which satishes fllll) . Dehne a 
new candidate for a Lyapunov function: 

(12) Vx { x ) xgM+. 

This function satishes Ra( 0) = 0) because (^'(0) = 0. Let us impose an additional assumption 
on the family {r'x)xew.+ of jump measures. 
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Assumption 3. There exists a Aq > 0 such that 

poo 

sup / < oo. 

x>0 Jo 

Under Assumption 3, we can dehne the following quantity for A G [0, Aq] and x G M+: 

1 

(13) K{x, A) = g{x)\ + -cr^(a;)A^ + j - l] u^{dy). 

Now comes one of the two main results of this paper. This hrst result is a statement of conver¬ 
gence with exponential rate, but it does not provide an estimate for this rate. 


Theorem 3.2. Under Assumptions 1, 2, 3, suppose there exists a A G (0, Aq) such that 
(14) lim A(a:, A) < 0. 

X^OO 

Then the reflected jump-diffusion is V\-uniformly ergodic, and the (unique) stationary distribu¬ 
tion 71 satisfies (vr, Vx) < oo. 


Proof. Note that Ux-uniform ergodicity and UA-uniform ergodicity are equivalent, because these 
functions are themselves equivalent in the following sense: there exist constants ci, C 2 such that 


V \ ix^ 

0 < Cl < , , < C 2 < oo for all x G M+. 

Vx{x) 

The corresponding reflected diffusion without jumps is totally irreducible, see [U Lemma 5.7]. 
The reflected jump-diffusion is also totally irreducible: With probability at least there are 
no jumps until time t > 0, and the reflected jump-diffusion behaves as a reflected diffusion 
without jumps. Therefore, by Proposition 13.11 it is sufficient to show that Ua is a Lyapunov 
function (in the sense of Definition [2]) for this reflected jump-diffusion. Apply the generator 
C = A + M from (El) to the function Ua. For x > S 2 , we have: V\{x) = V\{x). Now, the 
operator A from ([2]) is a differential operator, and its value at the point x depends on its value 
in an arbitrarily small neighborhood of x. Therefore, for x > S 2 , 


(15) 


AVx{x) = AVx{x) 


g{x)X + ^aflx)X^ 


Vx{x). 


Apply the operator J\f from ([7]) to the function Ua. For y G M+, we have: ip{y) < y, and 
therefore 

(16) Vx{y) = < Vx{y) = 

From (ITBl) . we get the following comparison: for x > S 2 and y G M+, 

(17) Vx{y) - Vx{x) < - e^^ = - l) = Vx{x) - l) . 

Because of ffT7|) . we have the following estimate for the operator M'. 

poo 

(18) AfVxix) < Vx{x) / - l] ufldy). 

Jo 

Combining ffTSj) and ffTSjl . and recalling the definition of k{x, A) from ffT3|l . we get: 

(19) CVx{x) = AVx{x) + AfVxix) < K{x,X)Vx(,x) for x > S 2 , Ag(0, Aq). 


Let us state separately the following technical lemma. 


Lemma 3.3. For every X G (0, Aq), the function CVx{x) is bounded with respect to x on [0, S 2 ]. 
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Assuming we already proved Lemma 13.31 let us complete the proof of Theorem 13.21 Denote 
(20) Co := sup [CVx{x) + koVx{x)] < oo. 

x£[0,S2] 

Combining (|20|) with ([19]), we get that 

k^Vx{x) < -koVxix) + Co1[o,s2](2^)) ^ ^ 
which completes the proof of Theorem 13.21 

Proof of Lemma The function Vx has continuous hrst and second derivatives, and by 
Assumption 1 the functions g and are also continuous. Therefore, AVx is bounded on [0, S 2 ]. 
Next, let us show that the following function is bounded on [0, S 2 ]: 


( 21 ) 


NV 




\y\{.y) - Vx{,x)\ = / Vx{,y)i^x{,^y) - Vx{,x)r{x). 


The function fl^ can be estimated from below by —Vx{x)r{x), which is continuous and therefore 
bounded on [0, S 2 ]. On the other hand, (l2l|) can be estimated from above by 


( 22 ) 


Vx(y)vMy) = e 


Xx 


oAy-x\ 


«Ax(d|/). 


From Assumption 3, it is easy to get that the function from fl2^ is bounded on [ 0 , 52 ]. This 
completes the proof of Lemma 13.31 and with it the proof of Theorem 13.21 □ 

Now, let us hnd some examples. Dehne the function 

pOO 

m{x):=g{x)+ {y-x)ux{dy), x e M+. 

Jo 

This is a “joint drift coefficient” at the point x G M+, which combines the actual drift coefficien 
g{x) and the average displacement y — x ioi the jump from x to y, where y ~ [r{x)]~^i'x{-)- One 
can assume that if m{x) < 0 for all or at least for large enough x G M+, then the process has a 
unique stationary distribution. This is actually true, with some qualihcations. 

Corollary 3.4. Under Assumptions 1, 2, 3, suppose is hounded on M_|_, and 

(23) lim m(x) < 0. 


Then there exists a A G (0, Aq) such that flTT|) holds. By Theorem 13.21 the reflected jump- 
diffusion Z = {Z{t),t > 0) is Vx-uniformly ergodic, and the stationary distribution vr satisfies 
( tt , 14) < 00. 

Proof. Because of Assumption 3, we can take the hrst and second derivative with respect to 
A G (0, Aq) inside the integral in flT3|) . Therefore, 

dK 


(24) 

(26) 

Letting A = 0 in 
(26) 


d\ 


= g{x) + (T‘^{x)\+ / P[y - x)vx{dy), 


d^K 

'W 

we have: 
dK 


= cr (x) + 


qKv P[y - xfvx{dy). 


d\ 


= g{x)+ {y - x)ux{dy) = m{x). 


A=0 
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By the condition ([23]), there exist mo > 0 and Xq > 0 such that 

(27) m(x) < —mo for x > Xq. 

There exists a constant Ci > 0 such that for all z > 0, we have: Applying this 

to z = \y — x\, we have: for A G [0, Ao/2], 

poo poo 

(28) / (2/-x)V("-"V,(d2/) <^1 / 

Jo Jo 

Combining fl28|) with Assumption 3 and the boundedness of a^, we get that the right-hand side 
of (|25|) is bounded for A G [0, Ao/2] and x G M+. Let C 2 be this bound: 

(29) C 2 := sup 

Ag(0,Ao/2] 


d^K 


dX^ 


x,X) 


By Taylor’s formula, for some A(x) G [0, A], we have: 
(30) 


BK X"^ r)^K 

K{x, A) = K{x, 0) + A—(x, 0) + y A(x)). 


Plugging fl26|) and K{x, 0) = 0 into fl30|) and using the estimate fl2^ . we have: 
(31) 

Combining (1^ with (I27|) . we get: 


K{x, A) < Am(x) -|- y A G 


0,^ 

’ 2 


C 

K{x, A) < K{X) := —Amo + y A^, x > xq, A G 

It is easy to see that K{X) <0 for A G (0, 2mo/C2]. Letting 

_ 2mo Ao 
'A — — 77;: — A —, 

C 2 2 ’ 

we complete the proof of Corollary 13.41 


0,^ 

’ 2 


□ 


Remark 1. In the setting of Theorem 13.21 or Corollary 13.41 the convergence of moments of 
P*(x, •) to the moments of 7r(-) follows from I/x-uniform ergodicity. Indeed, take an a > 0. 
There exists a constant C'(a, A) > 0 such that x" < C{a, A)V\(x) for x G M+. From (jHI) we get: 
for X G M+, t > 0, 


2 /“P*(x,d?/) 


l/“7r(d|/) 


< C{a,X)KVx{x)e 


—>ct 


4. Stochastic Ordering and Explicit Rates of Exponential Convergence 

In this section, we get to the main goal of this article: to explicitly estimate x, the rate of 
exponential convergence from (jH]). In case when the reflected jump-diffusion is stochastically 
ordered, and z = 0 in (|HD, we can (under some additional technical assumptions) get that x = fc. 
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4.1. General results for Markov processes. For two finite Borel measures v and v' on M_|_ 
with z/(M_|_) = we write v -< z/', or v' ^ z/, and say that v is stochastically dominated by 

z/', if for every z G M+, we have: z/([z, oo)) < z/'([z, oo)). 

Definition 5. A family {i'x)x>o of hnite Borel measures, with z/a;(R_|_) independent of x, is 
called stochastically ordered if z/ 3 ;(M+) does not depend on x, and z/j. ^ Uy for x < y. A 
Markov transition kernel P^{x,-), or, equivalently, the corresponding Markov process is called 
stochastically ordered, if for every t > 0, the family (P*(x, •))x>o is stochastically ordered. 

Remark 2. An equivalent dehnition of a Markov process X = {X{t),t > 0) on M_|_ being 
stochastically ordered is when we can couple two copies of this process starting from different 
initial points such that they can be compared pathwise. More precisely, for all x, y such that 
0 < X < y, we can hnd a probability space with two copies X^^'> and starting from 
= X and = y respectively, and X^^l(t) < X^^^t) a.s. for all t > 0; this follows 

from [T7] . 

In this section, we would also like to make 14 from (flU]) a Lyapunov function as in ([H]) with 
z = 0. However, we cannot directly apply the generator C from (|6]) to this function, for the 
reason we already mentioned: 14(0) 4 0, which contradicts ([3]). Neither can we use the function 
Vx from (fT^ : as follows from the proof of Theorem 13.21 we would have z = S 2 > 0, where S 2 is 
taken from ffTTi) . In [21], this obstacle is bypassed by switching to a (non-reflected) diffusion on 
the whole real line, but we resolve this difficulty in a slightly different way. The proofs of [2X1 
Theorem 2.1, Theorem 2.2], mainly use the Lyapunov condition (IH|) only “until the hitting time 
of 0”. To formalize this, let us adjust Dehnition HJ Let r(0) := inf{t > 0 | X(t) = 0}. 

Definition 6. A function V : M+ [1, cxo) is called a modified Lyapunov function with a 
Lyapunov constant A: > 0 if the following process is a supermartingale for every starting point 
A(0) = a; G M+: 

^iAr(O) 

(32) M{t) :=V{X{tAT{0))) + k H(X(s))ds, t > 0. 

Jo 

Remark 3. It is straightforward to prove that if 1/ is a Lyapunov function from Dehnition [2] 
with Lyapunov constant k and with z = 0 from ([HD, then 1/ is a modihed Lyapunov function 
with Lyapunov constant k. In other words, Dehnition [6] is a generalization of Dehnition [2] with 
z = 0. Indeed, because A4 is the generator of X, the following process is a local martingale: 

V{X{t))- [ A41/(X(s))ds, t>0. 

Jo 

If H is a Lyapunov function from Dehnition [2] with Lyapunov constant k and with z = 0, then 
the following process is a local supermartingale: 

M{t) := V{X{t)) - [ [-kV{X{s)) + 64o}(X(s))] ds, t > 0. 

Jo 

Moreover, this is an actual supermartingale, because it is bounded from below (use Fatou’s 
lemma). Therefore, the process {M{t A r(0)),t > 0) is also a supermartingale by the optional 
stopping theorem. It suffices to note that M(t A t(0)) = M(t). 

The following is an adjusted version of [2T1 Theorem 2.2], which states that for the case of a 
stochastically ordered Markov process with z = 0 in ([HD, we can take x = fc in ([HD- Note that 
we do not require condition (2.1) from [2T], but we require {j,V) < oo instead. For rehected 
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jump-diffusions, this assumption (tt, f/) < cxo can be obtained from Theorem 13.21 which does 
not state the exact rate k of exponential convergence. 

Theorem 4.1. Suppose X = {X{t),t > 0) is a stochastically ordered Markov process on R+. 
Assume there exists a nondecreasing modified Lyapunov function V with a Lyapunov constant 

k. 

(a) Then for every xi,X 2 G M+, we have: 

(33) ||P*(a;i, •) - P\x2, Oik < [V{xi) + V{x2)] t > 0; 

(b) For initial distributions /xi and /i 2 on M+, with {pi, 1/) < cx) and {p, 2 , V) < oo, we have: 

(34) IkiP* - /iaP'lk < [(hi, V) + (/i2, f")] t > 0; 

(c) If the process X has a stationary distribution tt which satisfies (tt, 1/) < oo, then this 
stationary distribution is unique, and the process X is V-uniformly ergodic with exponential 
rate of convergence k =k. More precisely, we have the following estimate: 

(35) \\P\x,-)-7i\\v<[{7i,V) + V{x)]e-'^\ t>0; 


Theorem 14.11 is an immediate corollary of Theorem 15.21 

4.2. Application to reflected jump-diffusions. To apply Theorem 14.11 to reflected jump- 
diffusions, let us And when a reflected jump-diffusion on M+ is stochastically ordered. 

Lemma 4.2. Assume the family (z^a;)xGR+ Is stochastically ordered. Then the reflected jump- 
diffusion from Section 2 is also stochastically ordered. 

Proof. This statement is well known; however, for the sake of completeness, let us present the 
proof. Let y > x > 0. Following Remark [2], let us construct two copies and of the 
reflected jump-diffusion, starting from Z^‘^l{f)) = X and Z^y\0) = y, such that a.s. for t > 0 we 
have: Z^^fit) < Z^^fit). Because z/ 2 ;(R+) = r{x) = r does not depend on a; G R+ (this follows 
from Deflnition [5]), we can assume the jumps of these two processes happen at the same times, 
and these jumps Ti < T 2 < ■ ■ ■ form a Poisson process on R_|_ with rate r. That is, — r„_i 
are i.i.d. Exp(r); for consistency of notation, we let tq := 0. Deflne these two processes on each 
[r„, Tn+i), using induction by n. 

Induction base: On the time interval [ro,ri), these are reflected diffusions without jumps on 
R+. We can construct them so that Z^^^t) < Z^^^t) for t G ko,ri), because the corresponding 
reflected diffusion on R+ without jumps is stochastically ordered. 

Induction step: If the processes are defined on [r„_i,r„) so that Z^^^t) < Z^y\f) for t G 
[r„_i,r„) a.s. Then by continuity Xn ■= Z^^\rn—) < yn ■= Z^yfirn—) a.s. Generate 2'0)(r„) ~ 
r~^u^^{-) and Z^yfirn) ~ so that Z^^\Tn) < Z^yfirn) a.s. This is possible by r'xni ) ^ 

iTy^f). Because the corresponding reflected diffusion without jumps is stochastically ordered, 
we can generate ZG) and Z^y'^ on (r„, r„+i) as reflected diffusions without jumps such that 
ZG)(t) < Z^y\t) for t G [Tn,Tn+i). Tliis completes the proof by induction. □ 

Next comes the central result of this paper: an explicit rate of exponential convergence for a 
reflected jump-diffusion on R_|_. 
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Theorem 4.3. Consider a reflected jump-diffusion Z = {Z{t),t > 0) on M_|_ with a stochastically 
ordered family of jump measures {yx)xm+- Under Assumptions 1, 2, 3, suppose for some A > 0, 

(36) -h'max(A) := supi^(a;, A) < 0. 

x>0 

Then the process Z is V\-uniformly ergodic with the exponential rate of convergence x = 
|-h'max(A)|. The (unigue) stationary distribution tt satisfies (vr, Vx) < oo. 

Remark 4. In certain cases, this exponential rate of convergence is exact: one cannot increase 
the value of x for the given norm for hxed A; see Lemma [6.11 in Section 6. 

Proof. That this process has a unique stationary distribution vr with (tt, Vaq) < oo follows from 
Theorem 13.21 In light of Lemma 14.21 to complete the proof of Theorem 14.31 let us show that 
Vao is a modihed Lyapunov function. Take an r; > 0 and let rijj) := inf{f >0\Z(t)< r]}. Let 
us show that the following process is a local supermartingale: 

(37) I4o(^(t A T{p))) + \K^M)\ / fAo(^(s))ds, t > 0. 

Jo 

Indeed, take a function Vx from flT^ with the function ip from flTT]) constructed so that S 2 < rj. 
From ([19]), we have: CV\{x) < K{x, A)KA(a^) fora; > rj. But K{x, A) < iFmax(A) < 0. Therefore, 

(38) CVx{x) < K^,flX)Vx{x) = -\K^.flX)\Vx{x), x > p. 

The following process is a local martingale: 

Vx{Z{t)) - [ CVx{Z{s))ds, t > 0. 

Jo 

By the optional stopping theorem, the following process is also a local martingale: 

_ rtr.T{ri) _ 

(39) VxiZit A T{p))) - / CVxiZ{s))ds, t > 0. 

Jo 

Observe that Vx{x) = Vx{x) for x > p, but Z{s) > p for s < rijf). Combining this with fl39D 
and fl38l) . we get that the process in fl37l) is a local supermartingale. It suffices to let r/1 0 and 
observe that t{p) f r(0). Therefore, for p = d the process in flTTll is also a local supermartingale. 
Actually, it is a true super mart ingale, because it is bounded from below (apply Fatou’s lemma). 
Apply Theorem 14.11 observe that the function Vx is nondecreasing, and complete the proof. □ 

The next corollary is analogous to Corollary 13.41 Its proof is also similar to that of Corol¬ 
lary 13.41 and is omitted. 

Corollary 4.4. Consider a reflected jump-diffusion on M+ from Section 2, with a stochastically 
ordered family (r'x)a:eM+ of jump measures. Under Assumptions 1, 2, 3, if 

snpm{x) < 0, sup cT^(a;) < oo, 

tc>0 a;>0 

then there exists a Aq > 0 such that iFmax(Ao) < 0, in the notation of fl56D . Therefore, 
the reflected jump-diffusion is Vx^-uniformly ergodic with exponential rate of convergence x = 
|A'max(Ao)|- The (unigue) stationary distribution'll satisfies (vr, Vaq) < cxo. 

Example 1. Consider the case when = 0: there are no jumps, this process is a reflected 
diffusion on the positive half-line. Assume 

sup 5 f(a;) = —g < 0, a{x) = 1. 

i>0 
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Then in the notation of fl36D . we can calculate 




supA) = sup g{x)X-\ - 

a:>0 x>0 2 



This function -ft'max(A) assumes its minimum value —g‘^/2 at A* = g. Therefore, this reflected 
diffusion is V^uniformly ergodic with exponential rate of convergence k = 'g^ j2. This includes 
the case of reflected Brownian motion on the half-line with negative drift from [211 Section 6]. 


Example 2. Consider the case g{x) = —2, a{x) = 1, and Ux = ^x+i for x > 0. In other words, 
this reflected jump-diffusion has constant negative drift —2, constant diffusion 1, and it jumps 
with rate 1; each jump is one unit to the right. Assumption 3 holds with any Aq. The negative 
drift “outweighs” the jumps in the positive direction: m{x) = —2 -|- 1 = —1. We have: 

A2 

K{x, A) = K{\) = -2X + — + e^ -1. 

For every A > 0 such that K{X) < 0, this process is lA-uniformly ergodic with exponential rate 
of convergence x = |A'(A)|. It is easy to calculate that K{X) < 0 for A G (0,0.849245). For 
example, the function K{X) attains minimum value —0.230503 at A* = 0.442954. Therefore, 
this reflected jump-diffusion is Va*- uniformly ergodic with exponential rate of convergence x := 
0.230503. Choosing larger values of A such that K{X) <0 (say, A = 0.8) results in lower rate of 
exponential convergence, but the norm which measures convergence becomes stronger. 


Example 3. Consider a reflected jump-diffusion with the same drift and diffusion coefficients 
as in Example 1, but with Vxiphi) = ^{y>x}^^~^(hl- other words, the jumps occur with rate 
z/ 2 .(M+) = 1, but each jump is to the right with the magnitude distributed as Exp(l). Then 
Assumption 3 holds with Aq = 1, and m{x) = —2 -|- 1 = —1. We have: 

K(x, \) = K(\) = -2A + A + ^ _ 1. 

This attains minimum value —0.135484 at A* = 0.245122. Therefore, this reflected jump- 
diffusion is VA,-uniformly ergodic with exponential rate of convergence x = 0.135484. 


5. The Case of Non-Stochastically Ordered Processes 

If the reflected jump-diffusion is not stochastically ordered, then we can still sometimes esti¬ 
mate the exponential rate of convergence. This is the case when this process is stochastically 
dominated by another reflected jump-diffusion, which, in turn, is stochastically ordered. 

Definition 7. Take two Markov processes X = {X{t),t > 0), X = {X{t),t > 0) with transition 
kernels {P*)t>o, (P) t>o on M+. We say that X is stochastically dominated by X if P*(x, •) ^ 
P (x, •) for alH > 0 and x G M+. In this case, we write X ^ X. 

The following auxiliary statement is proved similarly to Lemma 14.21 

Lemma 5.1. Take two reflected jump-diffusions Z and Z on M_|_ with common drift and diffusion 
coefficients g and which satisfy Assumption 1, and with families {i'x)xew+ o,nd {yx)xew+ of 
jump measures satisfying Assumption 2. Assume that Ux Fx for every x G M_|_, and the family 
(Fx)xgk+ is stochastically ordered. Then Z < Z. 

The next result is an improvement upon [211 Theorem 3.4]. We prove convergence in ll'Hv- 
norm, that is, uniform ergodicity, as opposed to convergence in the total variation norm, which 
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was done in [2T1 Theorem 3.4], In particular, if V = Vx, as is the case for reflected jump- 
diffusions, then we can estimate the convergence rate for moments, as in Remark [1] Such 
estimation is impossible when one has convergence only in the total variation norm. 

Theorem 5.2. Take a (possibly non-stochastically ordered) Markov process X = {X{t),t > 
0) which is stochastically dominated by another stochastically ordered Markov process X = 
{X(t),t > 0). Assume X has a modified nondecreasing Lyapunov function V with Lyapunov 
constant k. 

(a) Then for every Xi,X 2 G M+, we have: 

(40) \\P\x,, •) - Pfix 2 ,-fiW < [V{xfi + V{X 2 )] e-’^\ t > 0. 

(b) For initial distributions joi and /i 2 on M_|_ with (pi, R) < oo and (/i 2 , V) < oo, we have: 

(41) ll/iiP* - /i2Plu < [(/ii, R) + (/i2, R)] e-^\ t > 0. 

(c) If the process X has a stationary distribution tt which satisfies (tt, R) < oo, then this 
stationary distribution is unique, and the process X is V-uniformly ergodic with exponential 
rate of convergence k =k. More precisely, we have the following estimate: 

(42) ||P*(x,-)-vr||y < [(7r,R) + R(x)]e-'=*, t > 0; 


Proof. Let us show (a). We combine the proofs of pil Theorem 2.2], [221 Theorem 4.1], and 
modify them a bit. Without loss of generality, assume X 2 < Xi. Consider two copies Xi,X 2 of 
the process X, and two copies Xi, X 2 of the process X, starting from 

Xi( 0 ) = Wi( 0 ) = Xi, ^2(0) = ^2(0) = X 2 . 

Take a measurable function g : M+ —)■ M such that < R, and estimate from above the 
difference 

(43) |E^7(W(f))-E^7(X2(f))|. 

Because of stochastic ordering, using X 2 < xi, we can couple these processes so that 

(44) X2{t) < X2{t) < Xi{t), X2{t) < Xfit) < Xfit). 

Dehne the stopping time r(0) := inf{t > 0 | Wi(t) =0}. By (jHj), Xi(r(0)) = X2(r(0)) = 0, so 
r(0) is a (random) coupling time for Xi and X 2 : the laws of {Xi{t),t > r(0)) and {X 2 {t),t > 
r(0)) are the same. Therefore, 

E(7(Wi(t))lp>^(o)} = Eig{X2(t))l^t>T{o)}, 

and the quantity from fl43l) can be estimated from above by 

(45) |E5f(Xi(f))lp<^(o)} — E5f(X2(f))lp<:f(o)}I < E|5f(Xi(t))|lp<^(o)} + E|5f(X2(t))|l{t<T(o)}. 
Let us estimate the first term in the right-hand side of (I45|) . Because < R, we have: 

(46) E|5((Xi(t))|lp<:f(o)} < ER(Xi(t))lp<-(o)}. 

Next, because the function R is nondecreasing, 

(47) e'^‘ER(Xi(t))lp<^(o)} < e'=*ER(Xi(t))lp<^(o)} < E (Xfit A r(0)))] . 

Let us show that the following process is a supermartingale: 

(48) M{t) = [Xfit At)), t > 0. 


14 


ANDREY SARANTSEV 


Indeed, from fl32D . we already know that the following process is a snpermartingale: 

|•tAT{0) 

M{t) =V(Xi{tAT{0))) + k l/(Xi(s))ds, t>0. 

Jo 

Applying Ito’s formula to M{t) in fHS]) for t < r(0), we have: 

dM{t) = ke’^^V {Xi{t)) dt + e^Ml/ = e^*dM(f). 

This is also true for t > t( 0), because both M and M are constant on [r(0), cxd). Since M is a 
snpermartingale, it can be represented as M{t) = Mi(f) + M 2 (t), where Mi is a local martingale, 
and M 2 is a nonincreasing process. Therefore, 

M(t) = [ e^MMi(s) + [ e^MM2(s) =: Mi{t) + M2{t) 

Jo Jo 

is also a sum of a local martingale and a nonincreasing process. Thus, it is a local supermartin¬ 
gale. Actually, it is a true snpermartingale, because it is nonnegative (use Patou’s lemma). 
Therefore, 

(49) EM{t) < EM(0) = V{xi). 

Comparing fl46|l . flT7)l . fITH)) . fl49|) . we have: E| 5 f(Xi(f))|lp<:j=(o)} < V{xi). Similarly estimate the 
second term in the right-hand side of fl4^ . and combine this with (14^ : 

(50) \^9{Xi{t)) - Eg{X2m < [E(xi) -f ^( 0 : 2 )] e""*, t > 0. 

Taking the supremum in fl50|) over \g\ < V, we complete the proof of (I40|) . 

(b) Integrate over {xi,X 2 ) ~ /ii x /i 2 in fl5(l and take the supremum over \g\ <V. 

(c) Apply (b) to /ii = TT and /12 = Since V{x) > 1, we can take D = 1 + {ti,V) in ([9]). □ 

Now, we apply Theorem 15.21 to reflected jump-diffusions. 

Corollary 5.3. Take drift and diffusion coefficients g, cr^, satisfying Assumption 1. Take two 
families (^'x)a:GR+ (Px)xgk+ of jump measures which satisfy Assumptions 2 and 3, such that 
Vx Tx for every x G M+, and the family (Tx)x<=m.+ is stochastically ordered. Consider a reflected 
jump-diffusion on M+ with drift and diff^usion coefficients g^a"^, and the family {ox)x>o of jump 
measures. Let 

\2 roo 

K{x,\) = g{x)\ +a'^ix)^ + j - l]Vx{dy). 

Assume there exists a A* > 0 such that 

(51) supiC(a;, A*) =: iPmax(A*) < 0. 

x>0 

Then Z is V\^-uniformly ergodic with exponential rate of convergence k = |iPmax(A*)|. 

Proof. For each x G M+, the function y 1 —)■ — 1 is nondecreasing, and h>x ATx- Therefore, 

flOO flOO 

/ - 1] Uxidy) < / - l] Uxidy). 

Jo Jo 

This, in turn, implies that for x G M+, 

\2 roo 

(52) K{x,X^) = g{x)X^-^a'^{x)-;^-^ - l] z/a,(d|/) < iP(a;, A*). 

2 Jo 
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Comparing (1^ with fl52|) . we get: 

sup K{x, A*) < sup A*) < 0. 

3;>0 a:>0 

By Theorem l4.3l the process Z is VA^-uniformly ergodic, and its stationary distribution n satishes 
(^5 hA J < oo. Consider another reflected jump-diffusion Z with the same drift and diffusion 
coefficients g, cr^, and the family (■Fa;) 3 ;eK+ of jump measures. By Lemma ISTTl Z < Z. Similarly to 
Theorem l4.lL we can show that Vx, is a modified Lyapunov function for Z. Applying Theorem l5.2l 
and using Vx^ as a modified Lyapunov function, we complete the proof. □ 

Example 4. Take a continuous function -0 : M+ —)■ R+ such that 'ip{x) < x -1- 1. Consider a 
reflected jump-diffusion on R+ with g[x) = —2, cr^(x) = 1, and the family of jump measures 
{i'x)x> 0 i with Vx := 5^{x)- This family is not necessarily stochastically ordered (because -0 is not 
necessarily nondecreasing). However, ^ 6x+i, and we can apply Corollary 15.31 We have: 

_ \2 

iL(x, A) = —2A H —-—\- e^ — 1 

has minimum value —0.230503 at A* = 0.442954. Therefore, this reflected jump-diffusion is 
VA*-iRiiformly ergodic with exponential rate of convergence x := 0.230503. 

6 . The Best Exponential Rate of Convergence 

Consider a reflected jump-diffusion Z = {Z{t),t > 0) with constant drift and diffusion 
coefficients: g{x) = g, a^(x) = and with family of jump measures (iyx)x>o defined by 
I'xiE) = fi{{E — x) n R+) for E C R^, where /i is a finite Borel measure supported on R+. 
In other words, every Vx is the push-forward of the measure /i with respect to the mapping 
y ^ X + y. This process is a reflected Brownian motion on R+ with jumps, which are directed 
only to the right, with the magnitude and the intensity independent of x. 

Recall that the intensity of jumps originating from a point x G R+ is equal to r(x) := z/a;(R+), 
and its magnitude is distributed as \y — x|, where y ~ r“^(x)z/a,(-). In this case, the intensity 
of jumps is equal to r = /i(R+), and the magnitude is distributed according to the normalized 
measure This was the case in Examples 2 and 3 from Section 4. 

Then the family of jump measures {i'x)x>o is stochastically ordered. Next, 

K{x, A) = iC(A) =g\ + ^A^ + / - l] y{dz). 

^ Jm.+ 

From Theorem 14.31 we know that if 

(53) ^ / ^h(d2:) < 0, 

Jr+ 

then there exists a A > 0 such that K{X) < 0, and the reflected jump-diffusion is VA-uniformly 
ergodic with x = |iC(A)|. Actually, this rate of convergence is exact: one cannot improve this 
result. More precisely, for this A, one cannot And a x > |R(A)| such that the reflected jump- 
diffusion is VA-uniformly ergodic with exponential rate of convergence x. This is a counterpart 
of [21], Theorem 2.3], which finds the exact exponential rate of convergence in the total variation 
metric. Unfortunately, we cannot apply their results, because they require vr({0}) > 0 for a 
stationary distribution vr, which is not true in our case. As mentioned in Example 2, we can 
make a trade-off between the rate of convergence and the strength of the norm ll-Hv;,,. 
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Lemma 6.1. Under the condition fl53D . for every A G (0, Aq) such that K{X) < 0, and every 
Xi,X 2 G ffi+, we have: 

E,,VA(Z(t)) - E,,Vx{Z{t)) = {V^{x,) - VA(a;2))e-l'"W^ t > 0. 

Proof. Let x = |iL(A)|. We must go over the proofs of Theorems 13.21 I^TTl and 15.21 Using the 
notation from these theorems with an adjustment: fi{0) = we get: 

E,^Vx{Z{t)) - E,,V^{Z{t)) = EVx{X^{t)) - EV^iX^it)) 

= EVA(Xi(f))l{T-(o)>i} — EVA(W2(f))l{T-(o)>t} 

= EVA(Xi(t A r(0)))l{^(o)>t} - EVA(X 2 (t A r(0)))l{^(o)>t}- 
Multiplying by e^, we get: 

[E,,V^{Z{t)) -EMZm 

= E A r(0)))l{.(o)>i}] - E A r(0)))l{.(o)>o] 

= E A r(0)))] - E [e"(*^"®VA(X 2 (t A r(0)))] , 

because on the event {t > r(0)} we have: Xi{t A r(0)) = X 2 {t A r(0)) = 0. Next, if we show 
that 

(54) M{t) = A r(0))), t > 0, 

is a martingale for every initial condition Z’(O) = a; G M+, then the rest of the proof is trivial: 
just use E^M{t) = M(0) = Va(x) for x = xi and x = X 2 . Let us show that (|54|l is a martingale. 
We follow the proof of Theorems 13.2114.11 If a; > S 2 , where S 2 is taken from flTTD . then we have 
equality in f[T7|l . in flTK]) . and in flT^ . Indeed, take an a; > S 2 , and let y ~ r~^{x)h'x:{-)- Then 
1 / — X ~ p, therefore y — x > 0, and y > S 2 , (p{y) = y. Therefore, as in Theorem 14.11 the process 

rtAT{ri) 

Vx{Z{t A r(p))) + X / UA(Z(s))ds, t > 0, 

Jo 

is a local martingale for every 77 > 0, and hence for 77 = 0, because r(? 7 ) f t( 0). Similarly to 
the proof of Theorem 15.21 we can show that fl5T)) is a local martingale. Actually, it is a true 
martingale. Indeed, take an e := Aq/A — 1 > 0. Then for all x G M+ and f > 0, 

(55) E, sup [Ua(Z(s))]'+" < cx). 

0<s<t 

Indeed, we can represent Z{s) = B{s) + fi, where B = {B{s),s > 0) is a reflected 

Brownian motion on R+ with drift and diffusion coefficients g and cx^, starting from 5(0) = x, 
random variables fi ~ r“^/i(-) are i.i.d., J = (J7’(s),s > 0) is a Poisson process on R+ with 
constant intensity r, and B,J',^i are independent. Then 

J{t) 

(56) sup [Ua(Z(s))]^+^ = exp^Ao max 5(s)) exp^Ao 

0<s<t V 0<s<t j V ^ J 

— — 1=1 

The moment generating function G^{u) := Ee“^ of fi is hnite for u = Aq. Therefore, the moment 

generating function of the random sum of random variables is equal to 

J(t) 

G{u) := E exp ^77 = exp {r{G^{u) — 1)). 
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This quantity is also finite for u = Aq. Finally, Eexp (Aq niaxo<s<t-B(s)) < oo. Apply fl56|) and 
complete the proof of (l5^ together with the martingale property of fl5T|) and Lemma 16.11 □ 

7. Systems of Two Competing Levy Particles 

7.1. Motivation and historical review. Finite systems of rank-based competing Brownian 
particles on the real line, with drift and diffusion coefficients depending on the current rank 
of the particle relative to other particles, were introduced in [5] as a model for mathematical 
finance. Similar systems with Levy processes instead of Brownian motions were introduced in 
[35] . One important question is stability, do the particles move together as t —)■ oo7 We can 
restate this question in a different way: Consider the gaps between consecutive ranked particles; 
do they converge to some stationary distribution as t —)■ oo, and if yes, how fast? For competing 
Brownian particles, this question was resolved in jH Proposition 2], jSS] Proposition 2.2], [321 
Proposition 4.1]: necessary and sufficient conditions were found for stability. If the system is 
indeed stable, then the gap process has a unique stationary distribution tt, and it converges to 
TT exponentially fast. 

However, it is a difficult question to explicitly estimate the rate of this exponential conver¬ 
gence; it was done in a particular case of unit diffusion coefficients in [TT]. For competing Levy 
particles, partial results for convergence were obtained in [35] . However, an explicit estimate of 
the rate of exponential convergence remains unknown. In this section, we consider systems of 
two competing Levy particles. We improve the convergence conditions of [33] • In some cases, 
we are able to find an explicit rate of convergence. 

7.2. Definition and construction. Take a drift vector and a positive definite symmetric 
matrix 

(57) {g+,g.)eR^ A = 

Take a finite Borel measure A on Consider a Levy process L{t) = (L+(f), L_(t))', t > 0, on 
the space with drift vector {g^,g_y, covariance matrix A, and jump measure A. Take two 
real-valued r.c.1.1. (right-continuous with left limits) processes Xi{t), X 2 {t),t > 0, which satisfy 
the following system of equations: 

dXi(t) = 1 (Xi(t) > X2{t)) dL+(t) + 1 (Xi(t) < X2{t)) dL_(t); 
dX2{t) = 1 (Xi(t) < X2{t)) dL+{t) + 1 (Xi(t) > X2{t)) dL_(t). 

At each time t > 0, we rank the particles Xi{t) and X 2 {t): 

Y+{t) = Xi{t) V X2{t), Y_{t) = Xi(t) A X2{t), t > 0. 

In case of a tie: Xi{t) = X 2 {t), we assign to X 2 {t) the higher rank. (We say that ties are 
resolved in lexicographic order.) At this time t, the lower-ranked particle behaves as the process 
L_, and the higher-ranked particle behaves as the process L_|_. 

Remark 5. Assume there exist finite Borel measures z/_, z/+ on M such that 

A(dx+,dx_) = 5 o(cIt-i-) x z/_(da;_) -|- z/+(da;+) x 5o(da;_), 

then the jumps of the two particles occur independently. The jumps of the lower-ranked particle 
y_(t) and the higher-ranked particle H+(t) are governed by measures z/_ and respectively. 



a_|_ a _ 
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Lemma 7.1. There exists in the weak sense a unique in law solution to the system (1581) . More¬ 
over, consider the gap process Z(t) = l + (t) —Y_(t), t >0. This is a reflected jump-diffusion on 
M+ with drift and diffusion coefficients 

(59) g = g+— g- and cr^ = a+++ a_— 2a+_, 

and with the family {i'z)z&.+ of jump measures, where for every z G M+, the measure Vz is 
defined as the push-forward of the measure A under the mapping Fz : (a;+, X-) i—)■ |x+ — X- + z\. 


Remark 6. Here, we construct a slightly more general version of a system of competing Levy 
particles than in [35]. Indeed, in [35] they assume that the jumps are independent, as in 
Remark [5l moreover, the diffusion parts are also uncorrelated, a+_ = 0, and = i/_. 

Proof. Instead of writing all the formal details, which can be easily adapted from [5l |35], let us 
informally explain where the parameters from fl5^ come from and why the measure Vz for each 
2 ; G M+ is as described. Between jumps, the gap process moves as R+ — B_, where (i?_, R+) is a 
Brownian motion with drift vector and covariance matrix as in fl57|) . It is easy to calculate that 
R+ — B_ is a one-dimensional Brownian motion with drift and diffusion coefficients from (15^ . 

Next, let us show the statement about the family of jump measures. Assume that r is the 
moment of a jump, and immediately before the jump, the gap between two particles was equal 
to z: Z{t—) = z. Then i/+ = y+(r—) and y- = satisfy i/+ — i/_ = 2 :. Assume without 

loss of generality that Xi(r—) = y- and X 2 (r—) = i/+. (This choice is voluntary when there is 
no tie: ?/_ y+) but required if there is a tie: y- = y+i because ties are resolved in lexicographic 

order.) The displacement (x_,a;+) during the jump is distributed according to the normalized 
measure A, or, more precisely, [A(M^)] ^ A(da:+, da:_). After the jump, the positions of the 
particles are: 

Xi{T)=x_+y_, X2{T)=x+ + y+. 

The new value of the gap process is given by 

Z{t) = \X 2 {t) - Xi{t)\ = \x+ + y+- X-- yfl = \z + x+- x_| = Fz{x-,x+). 

Therefore, the destination of the jump of Z from the position z is distributed according to the 
probability measure Vz, which is the push-forward of the normalized measure 

[A(M2)]”^A(da:+,da:_) 

with respect to the mapping Fz. The intensity of the jumps of Z is constant and is always equal 
to the intensity of the jumps of the two-dimensional process L, that is, to A(M^). Therefore, the 
jump measure for Z{t) = z is equal to the product of the intensity of jumps, which is A(M^), 
and the probability measure z/^. The rest of the proof is trivial. □ 

7.3. Uniform ergodicity of the gap process. Lemma ITT] allows us to apply previous results 
of this paper to this gap process. First, we apply Corollary 13.41 Recall the dehnition of the 
function Vx from (fTOj) . 

Theorem 7.2. Assume that 

(60) ff e'^°^''^+'’''''^"'^A(da;+, da;_) < 00 for some Aq > 0; 


9+- 9 - + 



— X-] A(dx+, dx-) < 0. 


(61) 
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Then for some A > 0, the gap proeess is Vx-uniformly ergodic, and the stationary distribution tt 
satisfies (vr, Vx) < oo. 

We can rewrite the condition fl6T]) as < m_, where the magnitndes 


m+ = g+ + 


a:+A(da:+, da:_) and m_ = g^ + 


a:_A(da:+, da:_) 


can be viewed as effective drifts of the npper- and the lower-ranked particles: the snm of the 
trne drift coefficient and the mean valne of the displacement dnring the jnmp, mnltiplied by the 
intensity of jnmps. This is analogons to the stability condition for a system of two competing 
Brownian particles from |1]: the drift (in this case, the trne drift) of the bottom particle mnst 
be strictly greater than the drift of the top particle. 


Proof. Let ns check conditions of Corollary 13.41 Assnmption 1 and the bonndedness of are 
trivial. Assnmption 2 follows from the fact that the fnnction Fz{x+,x-) is continnous in z. 
Assnmption 3 follows from the condition flBU]) . Indeed, for (a;+,a;_) G and z G M_|_, we have: 
||a:+ — x_ + z\ — z\ < |a:_| -f |a;+|. Therefore, we get: for every z G M+, 

[ = [[ e^°ll"+-^-+"l-"lA(da:+,dx_) < //' e^°(l"-l+l"+l)A(da:+, da:_) < cx). 


Finally, let ns check the condition fl25]l . Indeed, 

/ [w — z]uz{dw) = // [l^;-|-— a;_| — ; 2 ] A(da;+, da;_) 


[a;_ — — 2^] A(da;+, da;_)-|- JJ [a;+— a;_] A(da;+, da;_) 

{x——xj^>z} {x——x^<z} 

[2a;_ — 2a;+— 2^;] A(da;+, da;_) + // [a;+— a;_] A(da;+, da;_). 

J Jr^ 

{x-—x+>z} 

However, 

JJ [2a;_ — 2a;+— 2^] A(da;+, da;_) < 2 JJ [x_ — x+] A(da;+, da;_). 

{x——x^>z} {x——x^>z} 

There exists a constant Cq > 0 snch that s < for s > 1. Therefore, for z > 1, 

[x_ — x+] A(da:+, da:_) < Cq JJ e'^°^^““'*+^A(dx+, dx_) 

{a;_—x+> 2 } {x-—x+>z} 

<Co JJ e^°(l^-l+l^+l)A(dx+, dx_) ^ 0 

{x-—x+>z} 

as 2 ; —)■ oo, becanse of flHHll . Combining this with previons estimates, we get: 


(62) lim / [w — zjuzidw) < // [x+— x_] A(dx+, da;_). 

Jr+ J Jr^ 

Combining (l6^ with (l6T]) . we complete the proof of 


□ 
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Remark 7. For the case of independent jumps, as in Remark [5], condition (l60|) can be written as 

/ oo roo 

e'^“l'^"lz/_(dx_) < cxo and / e'^°l'^+lz/+(dx+) < oo for some Aq > 0. 

-OO j —OO 

and condition flbTl) can be written as 

/ OO nOQ 

x+z/+(dx+) <g_+ x_u_{dx_). 

-OO J —OO 

7.4. Explicit rate of exponential convergence. In some cases, we are able to find an explicit 
rate x of exponential convergence for the gap process. This is true when the gap process is either 
stochastically ordered or dominated by some stochastically ordered reflected jump-diffusion. 
First, consider the case when the gap process is itself stochastically ordered. In addition to the 
assumptions of Theorem 17.21 let us impose the following assumption. 

Assumption 4. The measure A is supported on the subset {(a:+,a;_) G | > a:_}. 


Then the family of jump measures is stochastically ordered, because for x+ > X-, we 

have: \z + Xj^ — x_\ = z + x^ — x_, and this quantity is increasing with respect to 2 ;. Therefore, 
we can apply Theorem 14.31 We have: Fz{x+, X-) — z = \z + Xj^ — X- \ — z = Xj^ — X-. Applying 
the push-forward to the measure A, for g and from (15^ . we get: 

K{z, X)=g\ + - l] u,{dw) 


= g\+^\^ + 


^gA(x+-x_) _ A(da;+,da:_). 


This quantity does not actually depend on z, so we can denote it by K{\). If fc(A) < 0 for some 
A > 0, then the gap process is VA-uniformly ergodic with exponent of ergodicity x = |iF(A)|. 
In particular, if we wish to maximize the rate x of convergence, we need to minimize K{X). 
Actually, under Assumption 4, we are in the setting of Section 6, and x = |iF(A)| gives the 
exact rate of exponential convergence in the Vx-norm. 


Remark 8. For the case of independent jumps from Remark |5l Assumption 4 is equivalent to 
the condition that the measure is supported on M+, the measure V- is supported on M_, and 

2 poo poo 

K{X) = gX + yA^ + / [e^"+ - l] z/+(dx+) - 1 - / - l] u^dx.). 

Example 5. Consider a system of two competing Levy particles with parameters 

g+ = 0, g_ = 3, a++ = a__ = 1, a+_ = 0, 

and measures on R+ and v_ on M_ with densities 

z/+(da;+) = l{z:+>o}e~'^+dx+, iz_{dx_) = l{ 3 ,_<o}e^-da;_. 

In other words, the upper particle can jump upwards, and the lower particle can jump down¬ 
wards. For each particle, its jumps occur with intensity 1, and the size of each jump is distributed 
as Exp(l). Then conditions of Theorem 17.21 as well as Assumption 4, are fulhlled. Knowing the 
moment generating function of the exponential distribution, we can calculate 

K{X) = -3X + X^ + -^. 

1 — A 

This function obtains minimal value —0.0748337 at A* = 0.141906. Therefore, the gap process 
is Va*- uniformly ergodic with exponential rate of convergence x = 0.0748337. 
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The next example is when the gap process is not stochastically ordered, but is dominated by 
a stochastically ordered uniformly ergodic reflected jump-diffusion; we use Corollary 15.31 

Example 6. Take a system of two competing Levy particles with independent jumps, governed 
by measures z/+ = 0 and z/_ = 5i, with drift and diffusion coefficients 

9+ = 0, g- = 2, a++ = a__ = 1, a+_ = 0. 

As follows from the results of Section 6, the gap process is a reflected jump-diffusion with 
g = —2, = 2, and for x G M+. But ^ := S^+i for x G M+, and so 

K{x, A) = K{X) = -2A + A^ + - 1. 

This function assumes its minimal value —0.160516 at A* = 0.314923. Therefore, the gap process 
is Va*-L mifornily ergodic with exponential rate of convergence x = 0.160516. 
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